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Abstract 

The construction of needlet-type wavelets on sections of the spin line bundles over the sphere 
has been recently addressed in Geller and Marinucci (2008), and Geller et al. (2008,2009). Here we 
focus on an alternative proposal for needlets on this spin line bundle, in which needlet coefHcients 
arise from the usual, rather than the spin, spherical harmonics, as in the previous constructions. We 
label this system mixed needlets and investigate in full their properties, including localization, the 
exact tight frame characterization, reconstruction formula, decomposition of functional spaces, and 
asymptotic uncorrelation in the stochastic case. We outline astrophysical applications. 
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1 Introduction and motivations 

A growing literature has been recently concerned with (random or deterministic) spin functions, i.e. sec- 
tions of spin fiber bundles over the sphere (we defer a rigorous definition and more formal discussion to 
the next section). Actually the interest of the physical literature on such objects goes back for several 
decades, the seminal contributions going back to [32] and [ID]- In these papers, spin spherical harmonics 
were introduced in the language of physicists and used for the analysis of gravitational radiation. Much 
more recently, spin functions have found a crucial role in the analysis of cosmological observations, in 
particular in connection with so-called Cosmic Microwave Background polarization data. Polarization 
is a peculiar imprint characterizing the electromagnetic radiation which was emitted at the age of re- 
combination, some 13.7 billion years ago and in the immediate adjacency of the Big Bang; as such, it 
delivers information on a number of extremely important topics in the current landscape of physical and 
cosmological research, for instance on the existence of primordial gravitation waves, on the reionization 
era and on primordial non-Gaussianity. The literature on these issues is vast; we mention [3 IHl I48[ 124] 
for an introduction, while a massive amount of observations are currently being collected by satellite 
experiments by NASA and ESA ( WMAP and Planck, respectively). Spin fiber bundles will definitely be 
of the greatest interest for other areas of physical research in the next decade, for instance in the analysis 
of gravitational weak lensing on the images of galaxies [B] . We expect random sections of spherical fiber 
bundles to enjoy a growing relevance also outside the physical sciences, for instance in medical imaging. 
Despite such a rich environment from the physical sciences, the interest in the mathematical literature on 
spin bundles on the sphere has grown only very recently. In particular, some efforts have been entertained 
to extend to sections of spin fiber bundles the construction of spherical wavelets of a needlet type. Scalar 
needlets were introduced for the sphere by |40j , |41) ; the general case of compact Riemannian manifolds 
has been presented by [TH [TH |T7] . The analysis of the asymptotic properties of scalar needlets in random 
circumstances has been started by [2] ,[3], see also [3D], [T], [35], [S], [3], [IS] for further developments 
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and [33] , [55] , [H] , [II],[H],[3S],[3Z],[1I], among others, for applications to cosmological data. Spin needlets 
were introduced by |14| . in that paper, locahzation and uncorrelation properties are also addressed. The 
general case where the wavelet system need not be compactly supported in harmonic space is discussed 
by |18| ; stochastic properties and related statistical procedures are investigated in |13| , while applications 
to a CMB framework are provided by |,12 . In 1 , it is proved that spin needlets actually make up a tight 
frame system, with the same cubature points as the scalar case, and characterizations of Besov spaces on 
spin fiber bundles are discussed. Further results on the stochastic foundations of spin random fields are 
provided by [33^ and [32_j. 

The purpose of this paper is to consider an alternative construction for wavelets on spin fiber bundles. 
In particular, we focus on the case where the resulting needlet coefficients are (complex-valued) scalars, 
rather than spin quantities as in the spin needlet proposal. We label this system mixed needlets and 
investigate in full their properties, including localization, the exact tight frame characterization, recon- 
struction formula, decomposition of functional spaces, and asymptotic uncorrelation in the stochastic 
case; we outline also astrophysical applications. Concerning the latter, we stress in particular that mixed 
needlets allow for possibilities that were ruled out for the pure spin construction, such as the estimation 
of cross power spectra between scalar and spin components. 

The plan of the paper is as follows: in Section 2 we review some background material on spin fiber 
bundles, and we discuss some equivalent definitions which have been provided in the literature. In 
Section 3 we explain the construction of spin needlets, while Section 4 is devoted to the investigation of 
mixed spin needlets. Section 5 discusses characterizations of Besov spaces, establishing the equivalence 
of spin and mixed spin needlets in this regard and investigating on the properties of functional spaces 
for underlying scalar functions. Section 6 is devoted to directions for future research, with particular 
reference to statistical applications. 

2 Spin functions 
2.1 Some Definitions 

We begin by summarizing some definitions and basic facts about spin functions. For more details, the 
reader may consult our article [11] . 

We let N be the north pole of the unit sphere S'^, namely (0,0,1), and we let S be the south pole, 
(0,0,-1). Let 

Ui ^S^\{N,S}. 

If i?e 50(3), we define 

Ur = RUi. 

On Ui we use standard spherical coordinates (t?, ip) (0 < < tt, — tt < f < t^), and analogously, on any 

Ur we use coordinates (i?_r, ^Pr) obtained by rotation of the coordinate system on [//. 

At each point p of Ur, we let Pr{p) be the unit tangent vector at p which is tangent to the circle ^r = 

constant, pointing in the direction of increasing ipj^. (This is well-defined since i?N, i?S ^ Ur.) We think 

of Pr{p) as being our "reference direction" at p, relative to the chart Ur. 

If p e Ur^ n Ur^ , we let 

ippR^R-^ be the oriented angle from the reference direction Pr^{p) to pj^^{p). 

(See [M] for a precise explanation of which is the oriented angle. For example, at S, the oriented angle 
from i ti j is 7r/2; at N, the oriented angle from j to i is 7r/2.) 

Since R is conformal, the angle i^pR^R^ would be clearly be the same if we had instead chosen Pr{p) 
to point in the direction of increasing -Or, for instance. Thus ippRj measures "the angle by which the 
tangent plane at p is rotated if the coordinates are rotated by i?" . 
Now say f2 C §^ is open. Let 
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C^{n) = {F= (Ffl,)fl.g5o(3) : allFfl G C^iURnil), and for all i?i , i?2 e ^0(3) andallp G t/ij^ nf/fl, nl7, 
where V = V-pfeflJ- 

(Heuristically, a physicist would think of Fr as i^/ "looked at after the coordinates have been rotated by 
i?" ; at p, it has been multiplied by e**'^ , which is how physicists think of spin quantities behaving after 
rotation.) Equivalently, C^(J7) consists of the smooth sections over il of the line bundle with transition 
functions e''"^p«2«i from Ur^ to Ur^. 

Similarly we can define i^(ri) (the Fr need to be in L^). There is a well-defined inner product on 
given by {F, G) = {Fr, Gr); clearly this definition is independent of choice of R. 

There is a unitary action of 50(3) on L^(§^), given by F — ?> F^, which is determined by the equation 

iF^)iip)^FR{Rp). (1) 

We think of as a "rotate" of F. Thus we have two important relations: if we "rotate coordinates" , 
we have 

FR,ip)^e^^^FR,ip), 

while if we "rotate F" , we have 

(F«),(p) = Fr{Rp). 

Physicists would say that spin s quantities need to be multiplied by a factor e**'^ when rotated. 

2.2 Twisted bundles 

Although we will not make use of it in this article, it may be useful to recast the previous discussion in 
a slightly different form, as suggested by |33]- View SO{2) as a closed subgroup of S'0(3), with elements 
k G 5*0(2). This an Abelian subgroup, and assuming k is parametrized by the Euler angle 7 G [0, 27r] the 
irreducible representations of SO{2) are well-known to be one-dimensional and given by Ws{k) : C — >■ C, 
Ws{k)x= exp{isj)x, where s G N. Let g G 50(3), and consider the action 

k : {50(3) X C} ^ {S^C} , k{g,x) = {gk,exp{is-f)x) . 

We denote by £s the quotient space of orbits of the above action; that is, two elements (gi, xi) and (32, X2) 
belong to the same equivalence class if there exist k G 50(2) such that (52, 2:2) — {gik,Ws{k)xi). For 
s — 0, this is clearly isomorphic to {S^, C} , i.e. the space of complex- valued functions on the sphere. For 
s 7^ 0, we obtain indeed the same spin fiber bundle we defined before {Eg, t^, S^}, by taking the projection 

7r:£:, ->S2 ,7:ig,x)^gSO{2) , 

where we denoted as usual gSO{2) the equivalence class {gk, k G 50(2)} ; for 5 G 50(3), this is isomor- 
phic to 50(3)/50(2) ~S2_ 

2.3 Spin Spherical Harmonics 

Next we need some facts about the spin spherical harmonics; again the reader may consult [14] for further 
details. 

Let / G F^(S^) the space of square- integrable functions on the sphere; it is a well-known fact that the 
following spectral representation holds, in the F^ sense: 

f{x) = y^aimYimix) , aim = / f{x)Yim{x)dx , 
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or more formally 



1=0 



where {"H;} are the linear spaces spanned by the standard spherical harmonics {Yi^ 
are certain eigenfunctions of the (positive) spherical Laplacian 

Ag^Yira ^ l{l + l)Yi^ . 



-I < m < 1} , which 



Explicit expressions for the Yi„i may be found in [14] . 

We next define the spin-raising operator 9 and the spin-lowering operator 3. 

3, S are maps which take smooth spin s functions to smooth spin s + 1 (resp. s — 1) functions, and 
which commute with the actions of the rotation group ([!]). On a smooth spin s function F, we have 
{dF)R = Q^rFr, (gF)fl = SsrFr, where 



dsRFR{^,ip) 
hRFR{'dR,lfR) 



(sin 79k) 
{smi3R) 



d i d 

ddR siwdRdipR 
d id 



B-Or, sin d R dipR 
The spin s spherical harmonics, defined for / > |s|, are then given by 



(sini^ij) Fr{i9r,ipr) 
{sin 'd rY FR{-dR,(pR) 



(2) 
(3) 



(/ + .) 
(l + s) 



(BrYim , for s > 



1/2 



{-Qy'Yi^ , for s < , 



so that, if / > Isl 



[{i-s) {i + s + rf^^Yi^.s+i , 

-[{I + S) {I ~ S Yirn,s^i , 



(4) 
(5) 



see also [5T] and The llm.s, for I > \s\, —I < m < I, form an orthonormal basis for L^. In addition, 
Yim.,s is an eigenfunction of the (positive) spin spherical Laplacian 



A, = 



-m if s > 0, 



^33 ifs<0, 
acting on smooth spin functions, with eigenvalue 

eis = {l-\s\){l+\s\ + l). 



(6) 



(7) 



If s = 0, then As is just the usual (positive) spherical Laplacian. The formal adjoint of 3 (mapping 
smooth spin s functions on §^ to smooth spin s + 1 functions) is —3, so that A^ is formally self-adjoint 
on smooth spin s functions. 

For I > \s\, we let T-Lis denote the linear span of the Yim.s for —l<m<l. His is the eigenspace of A^ for 
the eigenvalue eis, and the direct sum of the His (for / > \s\, — / < to < /) is all of ^^(S^). 
We make some elementary observations about the eigenvalues e/^. 



If ; > |s|. 



eis ^ l{l + l) - \s\{\s\ + 1) < l{l + 1) ^ eio; 



(8) 
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and for any (always nonnegative) , 



(9) 



and for any / > and any s, 

ezo < fii+\s\,s- (10) 

Here ([5]) and (flUl) are trivial. To prove write e/+;'+i,o = e;o + e;/.o + 2(^ + + 1), then square both 
sides of (jni), to see that the inequality is equivalent to 



+ + 1) <{l + + 1), 

which is evident. 

Note that, for a spin s function Fg, we may speak unambiguously of the number [^^(a;)!, for any x e S^. 
We now prove the following inequality, by imitating the familiar method of proof in the case s = 0: 

Lemma 1 Say I > \s\, and that Y G His- Then 



\Y\\^<xl^^\\Yh. 



In particular, for any m, 



I -^m.s 1 1 oo — 



21 + 1 
47r 



Proof Let s+ = max(s, 0). Lr section 5 of [13] we showed that Zj^^, defined by 



zi,s = {-ly 



21 + 1 
An 



1/2 



Y, 



l, — s,s: 



(11) 



(12) 



(13) 



which we called the s-zonal harmonic for His, has quite similar properties to the usual zonal harmonic 
(for Hio)- Those properties, and ([T]) imply that the following argument is valid, just as in the case s = 0: 
Say p e §^ and choose R € 50(3) with RN = p. Then 



\Y{p)\ = \Y^iN)\ = \{Y^',Zi,s)\ < \\Y^4ZiA2 - \\Yh{Zi^s,Zi,sy^' = |1F|12|^^.(N)| V^. 



But 



\Zl,sim = I J2 YlrnA^)Yl^A^)\ 



21 + 1 



in fact, for any x €E S^, and any R' G 50(3), one has 



Ylm.sR'ix)Yirn,sR'ix) 



21+1 



(14) 



This completes the proof. 
For L>\s\, let 

L 

Vl.s = Hz. . 

l = \s\ 

We note the following Bernstein- type lemma, adapted from |19| . 



(15) 
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Lemma 2 A smooth spin s function F is in Vl ^ if o-n-d- only if there exist A > and B < e^+i^^ such 
that for every nonnegative integer N, 

\\iA,fFh<AB^. (16) 
Proof If F €E Vl,s, we surely have ^TE\\ . with A = H^'lb, B — eL,s, by the orthogonahty of the eigenspaces 

For the converse, say we have ([T6l) . Suppose I > L + 1, and Y E T-Li,s] it sufhces to show that {F^ Y) — 0. 
But, since is formany self-adjoint, for any N we have 

\{F,Y)\ = er/K^^, Af y)| = er/|(Af F,y)| < A{-f\\Yh. 

Since B < cl+i.s < e;^, this yields {F,Y) = upon letting N go to infinity. This completes the proof. 

From this we find the following important product property. (The case r ^ —s was first proved in 
[T] by developing the ideas of the subsection which follows. Here instead we adapt arguments from [19 .) 

Lemma 3 

VK,rVL,s ^ VK+L+\r+s\,r+s- 

Proof First note that the product of a smooth spin r function with a smooth spin s function is a smooth 
spin r + s function. 

Note next that if F^ G are respectively smooth spin p and spin q functions, then 

g(i^G) = {dF)G + F{m), (17) 

similarly for 3 in place of 9. p7|) follows at once from ([2]), once we note that for any p, R, as differential 
operators 

Spfl = Sofl +pcoi-dR, 

and similarly for g or p + g in place of p. 

To prove the lemma, it suffices to show that \i \r\ < k < K and |s| < ^ < L, and if F = Yfc^^.r and 
G = Yim,s for some /i,m, then FG G Vk+i+\r+s\,r+s- 

Iterating (|17p and the companion equation for 5, we find that we can write 

2N 

(-A.)^(FG)^^T„ (18) 

where 

each Tj is a sum of (2^) terms of the form {DF){D'G), where 
I? (resp. D') is a j-fold (resp. (2A^ — j)-fold) product of 3's and 9's, in some order. 

Note that the constants appearing on the right sides of ([2]) and ([3|) are equal to zt^/ei^ or ±y/e/~I7. Thus, 
by (H)), if Z? is as above, 

DF = or bYk^^^r' for some r', where |5| < e^Q . 
Thus, notation as above, by Lemma [I] we have 



\{DF){D'G)\\2<\\DF\\2\\D'G\\^< 



2? + 1 972^(27V-j)/ 



47r 



2 
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so that, by the binomial theorem, 



An 
where 



2' + 1 j^2N 



by dS]). Set B = Bl By ([TOl), B < efc+i+|r+s|+i,r+s- By Lemma[2l we see that FG G Vk+i+\r+s\,r+sj 
desired. 



as 



2.4 Connection with Wigner D matrices 

Next, we shall explain the connection of spin spherical harmonics with Wigner D matrices. This connec- 
tion provides an alternative point of view, but it is not necessary for the rest of the article. For further 
details on this connection, the reader may consult [13] and [33] . 

It is well-known that the elements DI^q, m — —I,..., I of Wigner 's D matrices are proportional to the 
standard spherical harmonics Yim ■ It turns out that this equivalence holds in much greater generality, in 
fact one has that (compare |33| for a discussion of phase conventions) 



2^ + 1- 



Yir^A^,^) = {-i)"^' ^=l—^Dl,^_,i^,^,0) . (19) 

Here in place of Yim,s, ^, we should have written Yim.sR, fR, throughout, but we drop the reference 
to the choice of chart for ease of notation whenever this can be done without the risk of confusion. 
Many of the properties of spin spherical harmonics follow easily from their proportionality to elements 
of Wigner's D matrices. Indeed, for instance, their orthonormality 

Yim^s{p)Yi'm',s{p)dp = / / Yi,n,s{-i},f)Yi'm'{^,'p)sin'dd'ddip^S\S'^ ; 
§2 Jo Jo 

is immediate. Also, viewing spin-spherical harmonics as functions on the group 50(3) (i.e. identifying 
p = (t9, lyj) as the corresponding rotation by means of Euler angles), and using and the group addition 
properties we obtain easily 

' 2^ + 1 



Yi^^s{p)Yi„,^,ip') = ——Y^Di„^_,iif,§,0)Dl^^,iip',^',0) 

m= — l m 

21 + 1 



An 



DiM{p,p')) , 



where Tp{p,p') denotes the composition of the two rotations (explicit formulae can be found in '49'). In 
the special case p — p\ we recover ([T4| . 



2.5 E and M modes 

For a smooth spin function on we have the expansion 

F = Y.Y1 aim.sYim^s ■ (20) 



I m 
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with rapid decay of the ,j in I. From (j20p . a further, extremely important characterization of spin 
functions was first introduced by [42^, see also [lOj and [14j for a more mathematically oriented treatment. 
In particular, it can be shown that there exist a scalar complex-valued function 



g{d,^)^n{g} + t^{g} , (21) 



such that. 



where 



'^^O.lm-EYlm,s + i'^^aim-lVlYlm^s 1 (22) 
Iva Im 



F^ = {^y^{g} ,F^ = {^y^{g} . 
Note that a/„i,s = aim-.E + iaim-M, where aiTn;E = ai,-m;E, aim-M = ai-m-M. It is also readily seen that 

0-1771,3 + W~^m^ = 0,i„,.E + iahji-M + 0-ijn;E ~ io,i,n-M — '^^ij^E , 
ai77i,s — ai-m,s — aim-.E + iai77i;M ~ aim-E + "iciini-M — 2iaim-]\I . 

In the cosmological literature, {aim-E} and {a;,„.M} are labelled the E and M modes (or the electric and 
magnetic components) of CMB polarization. 



3 Spin needlets 

We now recall the construction of spin needlets, see [14], [12], [18], [13] and [I] for further details and 
discussions. Fix a "dilation parameter" _B > 0; B is often chosen to be 2, but it is sometimes useful to 
let it take other values. Let cj) he a C°° function on R, symmetric and decreasing on supported in 
1^1 < 1, such that < (/)(C) < 1 and (f){^) = 1 if |^| < Let 

bHO^Hi)-HO>0 (23) 

Note that supp6 C [1/B, B], and that 

^52(^)^ lim 0(X)^lforanC>O. (24) 

J 

Of course the sum on the left side of (|24p is zero if ^ = 0. 

Let T be a positive self-adjoint operator on a Hilbert space H, and let P be the projection onto the null 
space of T. It is a special case of Theorem 2.1(b) of [15], that we may use the spectral theorem to replace 
f by T in ([HI), obtaining that 

J2bi^f = I-P. (25) 
j 

where the sum converges strongly. 

We take n = T = A„ Pn = H,,,,,. Thus, if = E; E™ aim,sYim,s G n, then 

^(^)^ = E E bi^)oi7nA,s. (26) 

/ m 

From this, it is easy to check (1^51) directly. Note also that PH is finite-dimensional (in fact, 2|s| -I- 1 
dimensional). Note moreover that b{^^) = for j sufficiently negative, specifically if B^^ < e|s|_|_i the 
smallest positive eigenvalue of As. 
For X e Uji, let 

A,(x, y,R)^J2Yl b{^)Yi7n,sR{^Wi^Ay) ■ (27) 

I m 
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Then evidently, if F{y) = Y.i Y.m "-im^sYimAy) ^ ^^^ve 

I m '' 

Here the integral is over It is important to note that, in the notation of subsection 12.31 

b{^)F e Fi^,^,, (29) 
where Ljs {'■= Lj if s is understood) is the largest integer with cl^s < i?^'-'^^'. In particular 

L, ^ (30) 

as j — >■ cxD. 

Now, take F ^ {I — P)H. Apply both sides of (^5]) to F, and take the inner product with F. We find 

\\F\\h - E \\bi^)nh = T.f \b{^)m^)dx , (si) 

while, as long as x e Ur, 

\K^)F?ix) - I / A,{x,y^R)Fiy)dy\' . (32) 



By ([^ and Lemma [31 \b{^)F\'^ £ ^^2Lj,o- That is, it is the restriction to the sphere of an ordinary 
polynomial of degree at most 2Lj ^ . Accordingly, as noted by [1] it is possible to follow the method 
used in [40] in the case s = 0: By familiar results for polynomials on the sphere, then, there is a constant 
c > 0, such that for each j, there is a c/ {Lj + l)-netEl jfc} of points on the sphere, and cubature weights 
{Ajfc} ^ {Lj + 1)~^, such that for every polynomial q of degree at most 2Lj, 

/ q{x)dx = ^ Xjkq{(.jk) (33) 
k 

Thus, for F £ {I — P)'H, and provided ^^j, £ Rjk, we in fact have 

W^Wh = EE^^-''-! / mjk,y,Rjk)Fiy)dy\^ . (34) 

j k 

In other words, for F £ {I — P)Ll, 

ll^lli?=EEl(^>.M)^ (35) 

j k 

where the spin needlets '4'jk,s ^'''^ defined by 

^jkAy) = V^^Mjk^y^R]k) = y^^K^^^)^Yi<nAy)Yini,sF(.,A£.jk) ■ (36) 

I ni 

Since ei^i ^ = 0, each ipjk,s ^ ^ P)^^- Consequently the {ipjk.s} '^''6 tight frame for (/ — P)L^- 
For F as above, we also define its spin needlet coefficients by 



l^jk^s = {F, ^jk,s) = V^E ^(~^) E a'm.silm.sfl,, i^jk) ■ (37) 



B3 



^See e.g. [3] for the definition of e-net. 
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By general frame theory, ii F E (I — P)Ll, we have the reconstruction formula 

^ = EE/^A>As- (38) 

j k 

Remarks 1. The choice of Rjk does not affect any of the terms on the right side of (|35|) or ((38)) . For 
this reason, and for simphcity we will sometimes omit the Rjk subscript in the formulas (j36l) and (I37p 
for ipjk,s ^^'^ Pjk,si when this can be done without causing confusion. 

2. We are ignoring the finite-dimensional space PLil- This is acceptable, because in astrophysical appli- 
cations, the interest is in high frequencies. 

3. One can use more general 6, than those we used in (P^. to construct spin wavelets on the sphere, as 
in [M], [18]. This leads to nearly tight frames with other interesting properties. For instance, one can 
arrange for the support of the frame elements at scale to be contained in a geodesic ball of radius 
CB~^ (for some fixed C). 

4. We will use the following notation and observations in Section 5. Let us set 

Q, = h{^f. (39) 
Following the arguments of ([5^ - ([55)) . but now without summing over j, we have that for F E L^, 

(Q,F,F)=^|(^^,V^^.,,,)P, (40) 

k 

After polarizing this identity, we see that for F E L^, 

g,.F = ^(F,V,,,.)V',M- (41) 

k 

In (iMl), b{^) = unless ^ E {l/B,B). Thus, for all j, 

i'jk^s^VL.eVL,.,. (42) 

Accordingly, 

Qj : Ll ^ Vl, e , and Q, = on [Vl, G J^- (43) 

For any integer N , let 

TV 

Pn= (44) 

j=-oo 

As we know, the sum in pi)) is actually finite, and, by (|25p . 1 ~ P strongly as — cx). By this, 

(|42|. and (gS]), we have 

^jk,s = (<3j-i + <3j + Qj+i)V'jfe,s, (45) 

for all j, fc. 

4 Mixed needlets 

We now present a construction of a different tight frame for spin functions, which we shall call mixed 
needlets. 

For now we work with ordinary functions. Let r be a fixed integer. Let be the space of spherical 
harmonics on §^ of degree Z, and let W = 0;^|r| '^i ^ We consider (1^ with 'K ~ W and with 
T replaced by T = A — |r|(|r| -I- 1), so that PW = 'H\r\- Since ew — e;o — + l)i '^'^ have that 

TYim = eirYim for Z > \r\. 
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In this situation, (^5]) leads to evident modifications of ([M)) . To modify the equations, we simply 

replace T by T, and take s = 0, the sole exceptions being that we write ew^ el instead of eio, bl qO- Of 
course we may disregard all the rotations R etc. We find that {rj^f.} is a tight frame for (/ — P)7i, where 



V.kiy) = E ^(^) E Ylm{y)y imi^jk)- (46) 

l=\r\ m 



The last step in constructing mixed needlets is to change the variable name r to s in (|46p . then to note 
that 'H'^ is unitarily equivalent to L^, by means of the unitary equivalence U, where U{Yim) = i^m.s (for 
I ^ \s\)- Thus, if P is now, once again, as in the previous section, so that PL^ — T-L^si^s^ then the mixed 
spin needlets (or mixed needlets for short) {"^j/t sm} ^ tight frame for (/ — P)Ll, where 



l=\s\ m 



For F — "^rn 0'im,sYim,s G L^, wc also define its spin needlet coefficients by 



Pjk,sM = '^3k,sM) = E ^(~^) E "■hn,sYhn{S.jk) ■ (48) 

By general frame theory, if F G {I — P)L'j., we have the reconstruction formula 

-P^ = E E l^3k,sM ^jk,sM ■ (49) 
j k 

In this "mixed" situation, we set Qjm — Ub{^^^)'^U~^ . Then "mixed" analogues of (|i(7)) - (H5|) hold; it 
is only necessary to replace by Qjm and V-'jfe.s by '4'jk,sM in those equations. 

But in fact we also have: 
Lemma 4 

QjM — Qj- 

Proof It suffices to show that these bounded operators agree on the orthonormal basis elements Yim,s- 
But 

QjMYlm,s = Ubi^fU-^Yi„,,s = Uh{^fYim - Ub{^fYun = K^fYl^^s = QjY^,. 

as desired. 

In brief, the construction of spin needlets proceeds by applying the methods of ([25)) - (p4|) to &(-^)^, while 
the construction of mixed needlets proceeds by applying the same methods to the unitarily equivalent 

operator fc(-^P)^, then invoking the unitary equivalence. Of course this unitary equivalence is only 
effective on L^, so there is no evident reason to think that there would be an effective theory in other 
function spaces for mixed needlets. However, the main point of this article is that mixed needlets do 
have nice mathematical properties beyond the theory, as well as useful astrophysical applications. In 
particular, as we shall show in the next section, they satisfy the usual needlet near-diagonal localization 
property, in the same sense as spin needlets were shown to in [T4^. 

To understand better the meaning of mixed needlets and its relationship with the existing literature, we 
start from (|22p . and introduce the notation 

fEix) C/-iF^ = E"'™^^^'™(^) ' /*^(^) •= U-'F^' = E°'™;^^^'™(^) • (^0) 

Im Im 
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Clearly and /m are well-defined scalar functions whicli are uniquely identified from Fs] as recalled 
above, in the s = 2 of interest for the physical literature they are labelled the electric and magnetic 
components of the spin field (in the physical literature, /m is rather written /b, but we already devoted 
the letter B for another purpose). Of course, it is possible to implement a standard (scalar) needlet con- 
struction on these spaces, enjoying the well-known properties of needlets (and indeed the same argument 
could be considered for other spherical wavelets). The interesting question to address is clearly what are 
the properties of such a procedure when viewed as acting on the original spin space Lg. 
More precisely, a direct idea to implement a wavelet transform on a spin random field would be as follows. 
Start by to evaluating the spin transforms 

where 

1 r _ T 1 r - 1 

aim.s = 2 {O'lm.s + 0-l,-m,sf + {aim,s — ai,-m,sf = aim;E + ^aim;M ■ 

Note however that it is not true that Re{aim:s) — o-im-.E, I'mio.hn-.s) = o,im;M, indeed aim-.E, o,im:M are 
complex-valued, and we have 

1 r _ T 1 r • - - 1 

aim-E = -^{aim^s + ai^m.si — -^{aim-E + iaim-^M + ai^^m^E — iai,-rn,M } , 

if _ T i r . _ .1 

aim;M — -~-^{ahn,s ~ ai-m,si = —-^{aim-E + ia,irri-M — ai^-m.E + iaim-M} , 

where we use the (involutive) property aim-E — ai,~m;E, aim-.M — aL-m-M- This property uniquely iden- 
tifies the spherical coefficients aim-E, ciim;M ■ Note that aim,s is involutive if and only if the M component 
is identically null, while if and only if the E component vanishes iaim.s is involutive. 
It is then readily seen that 



/3 



jk,M 



^' Im 



Im 



— Pjk;E + iPjk;Mi (51) 

where f3jk;Ei Pjk-M are real, i.e. (3^,,.^ = Re{/3jk,M)^ f^jk;M = I'm{Pjk,M)' because aim-E,aim-M and 
Yimi^jk) are involutive. It is immediate to verify that i3jk-^E^Pjk,M are exactly the coefficients we would 
obtain by evaluating a very slightly modified standard (scalar) needlet transform on the scalar functions 
Iei Im- (It is slightly modified because one would be using the rj^j. of (f46|) with r = s instead of the usual 
needlets, for which one would use h{l/B^) instead of h{^/eu/ B^) in 



4.1 Localization 

We fix an integer s. Because of (j36|) . localization properties of the spin needlets ipjks may be derived 
from localization properties of the kernel {x,y, R), defined in (j27l) . Often we will only need information 
about its absolute value |Aj(x,?/)| (which we write as shorthand for \Aj{x^y, R)\ for any R). 
For a full understanding, we need to consider the kernel 

Hx, y, t, R, R\ 5) = X! X! 9{tVei?)yim,sR{x)Y hn,sR' (y) , (52) 

/ m 

for t>0, R,R' e 5*0(3), g G (M). In particular 



\Aj{x,y)\ = \K{x,y,B-\R,R\h)l 



(53) 
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for any R, R' . 

When g is fixed and understood we will write A(a;, y, t, R, R') for A(a;, y, i?, i?', g). In the case s = 0, A 
does not depend on i?, i?', and we simply write A(a;, y, t) for A(x, y, t, R, R'). 

In the case s = 0, localization properties of a variant of this kernel (where e;o in (l52t is replaced by P) 
were derived in [4^, [41]. For the kernel as it stands in ([27| . as well as analogous kernels on smooth 
compact Riemannian manifolds, localization results (including results where b need not have compact 
support away from 0) were proved in [15 . 

For general s, the localization properties of A were proved in [T3] and |18) . 
For the purposes of this article, the relevant localization results are: 

Lemma 5 Say s = 0. 

(a) Suppose g G C^(0,oo). Then for every pair of C°° differential operators X {in x) and Y (in y) on 
S^, and for every integer r > 0, there exists C > as follows. Suppose degX = j and degy = k. Then 

\XYA{x, y,t)\< —-Yr, (54) 

|l + fe^} 

for all t > and all a;, y G S^. 

(b) If A > is fixed, and g £ C|?°(R) is even, then the conclusion of (a) remains true as long as < t < A. 
Remarks 1. In relation to the function / of [TS], our g{(^) = /(5^). 

2. Part (b) was shown in [T5] for < t < 1 - see the last paragraph of section 4 of that article. For 
t £ [1,j4], the estimate ((M)) is trivial, for then the right side is uniformly bounded below by a positive 
constant, and the left side is uniformly bounded above by a positive constant, as is apparent from (152^ .) 

3. Say [p,q] C (0,cx)). Using the remarks preceding Theorem 6.1 in [M], one sees that the constant C 
appearing in ()54|) may be taken to be uniform for g ranging over any bounded subset of the Frechet space 

{[p, q]). Indeed, this follows easily from Lemma [5] and the closed graph theorem. 

Lemma 6 Say s is a fixed integer. 

(a) fl4;l Suppose g £ C^(0,oo). Say R,R' £ SO{3). Then for every pair of compact sets J-ji C Ur and 
J^R' C Ur', every pair of C°° differential operators X {in x) on Ur and Y {in y) on Uri, and for every 
integer r > 0, there exists C > as follows. Suppose degX — j and degl^ = k. Then 

\XYA{x,y,t,R,R')\<- — -^r, (55) 

for all t > 0, all X £ Fr and all y £ Fri . 

(b) fT8 l If A > is fixed, and g £ C^(M) is even, then the conclusion of (a) remains true as long as 
<t< A. 

Remark Note that if X = Y = the identity operator, \XYA{x,y,t,R,R')\ is independent of R,R', and 
then ([55]) holds for aU x,y £ 

We are now going to show that a similar localization result to Lemma [S] (a) holds for the mixed needlets. 
Consider then the "mixed needlet kernel" 

AM{x,y,t,R) = b{t^/en)Yim,sR{x)Yi^{y) ; 

Irn 

It should be observed that 

Am {x, x,t) = for all x e S^ i > 0, (56) 

because 



b{tV^)^im,sR{x)Yim{a 



Im 



E Htv^s) d\^^j,{x)d'^^{x) 



= , 
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by the unitarity properties of the Wigner's D matrices. In contrast, in the unmixed situation |A(a;, x,t)\ = 
X); ^(^\/e;7)|, by (ITU) . If ^ 6 > 0, as is usuahy the case in appUcations, this is nonzero. Qual- 
itatively, then. Am is different from A; (|56l) might even lead one to suspect that Am might not be 
well-localized. However, the methods of our article [IT show that it is: 

Lemma 7 Say s is a fixed integer. 

(a) Suppose b £ C^{0,oo). Say R,R' G 50(3). Then for every pair of compact sets C and 
J'R' ^ Ub.', every pair of differential operators X (in x) on Ur and Y {in y) on S^, and for every 
integer r > 0, there exists C > as follows. Suppose degX — j and degY = k. Then 

\XYAM{x,y,t,R)\ < \\r , (57) 



2 



for all t > 0, all X £ Tr and all y G § 

Proof. We shall modify the proof of Theorem 6.1 of [T3]- To make that easier, we set /(^^) = b{Oj so 
that / S C^(0, oo); say that in fact supp/ C [p,q], where p > 0. We assume s > 0; similar arguments 
will apply when s < 0. Of course the case s = is handled by Lemma [5] 
We have 

Am {x, y, t,R) f{t'^eis)ICMB.{x, y) , 



1>S 



where 



Thus 



where 



and here 



/^Mni^^y) = ^Ylm^sR{x)Yi,n{y) 



Am {x, y,t,R) = ^ ^(^2^^^ ) ^ j IL^^ y) 



1>S 



(i + sy. 



l^o{x,y) = ^Yim{x)Yira{y) 



in the x variable. 

As in the proof of Theorem 6.1 of [13], we note that 



ji + sy. 

where 7^ := k{k - 1). As in [H] we choose To,Ti > with j^Tq < p/2, e^+i^sTj > q. We note that 
^M{x,y,t,R) = for i > Ti, and that ((57|) is trivial for t in the compact interval [Tq,Ti\ C (0,cx)). 
(Indeed, there the right side of ([57]) is uniformly bounded below by a positive constant, and the left side 
is uniformly bounded above by a positive constant.) It is then enough to focus on t G (0,To], We now 
define 

fiu-s{s + l)t^) 
Jt{u) := — , 

supported in the fixed compact interval [p, qi] := [p, q + s{s + l)!]^ ]; we note that the denominator does 
not vanish in the interval. Then, using ([5]) and (j58p . we write 




l>s V ^" ■ l>s 
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for 



1>S 



Now note that the functions ft for t S (0,To] form a bounded subset of C^([p, gi]), and recah Remark 
3, after LemmaO Choose smooth differential operators on in x, of degree s and j, which agrees with 
S^jj, and X respectively, in a neighborhood of Tr. As in the proof of Theorem 6.1 of [13], we find 



\XYAMix,y,t,R)\=t' Xrg^'^A[,](a;,y) 



{l + dix,y)/ty {l + d{x,y)/ty 



as desired. 



Remark 8 For astrophysical applications, it is very common to observe spin random fields only in a 
subset of the sphere, i.e. §^\G, say, where G C §^ is a region contaminated by foreground emission, for 
instance the Milky Way radiation. Lemma^ implies that, for t — 1,2, ... 



jk-E 



Re- 



Fs{x)ipjk,sMdx 



< 



< 



Ps{x)'il}jk,,j^dx 

Ps{x)\ l^jk^sAildx 



Fsix)ipik,sMdx 



< <isup|F,(2;) 



{l + Bidii^„G)y 



^jk;M 



— Im 



Fsix)'ip.k,sMdx 



< 



< 



Im { i/S.k-M 



Fs{x)')pjk,sMdx 



sup \Fs{x)\ 
xeG 



{l + Bid{i^^,G)y 



where denotes Lebesgue measure on In other words, for all the coefficients corresponding to 
locations ^jj, G §^\Ge (where G^ := ^ G §^ ; G) > e) the mixed needlet coefficients are asymptotically 
unaffected by the presence of unobserved regions. This is clearly a property of the greatest importance for 
cosmological applications. 

We now derive the following corollaries from Lemmas |6] and [71 which will be essential for the characteri- 
zations of Besov spaces in the next section. 

Corollary 9 (a) In (a) or (b) of Lemma\^ we have that for some G > 0, 

\K{x, y, t)\dx < G, / \A{x, y, t)\dy < C, 



where the integrals are over §^ . 

(b) In Lemma^ we have that for some G > 0, 



\AM{x,y,t)\dx < G, / \AM{x,y,t)\dy < C, 



where the integrals are over 
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Proof This follows at once from the fact that J[l + d{x,y)/t] ^ dx < Cpft" for any N > n (see for 
example, the third buUeted point after Proposition 3.1 of |15|). 

Corollary 10 Say 1 < p < oo. 

(a) For each j, Qj,Pj ■ — ?► Lf. If p > 2, the restriction of Qj,Pj to C is bounded on L^. If 
p <2, Qj,Pj may he extended from to be bounded operators on Lf. 
Further, the operators Qj,Pj are uniformly bounded on for —oo < j < oo. 

W ¥3k,s\\^^\^3k,sM\\^<C2-^■ 

Proof For (a), recah from ^ that Qj = b{^f on i^.Take A as in (l52|) for g = h'^. Then for x e Ur, 
y £ Uri , we have 

iQjF)Rix) - J A{x, y, t, R')Fr, {y)dy (59) 

so that 

\{Q,F){x)\ < J \A{x,y,t)\\F{y)\dy 

for F e L^. Part (a) for Qj is apparent from this and from Corollary IHl^a). (Note also that ([5^ continues 
to hold for all F G (one uses a density argument if p < 2).) Similarly, part (a) for Pj also follows from 

Corollary |9] (a) (where one references part (b) of Lemma|6]), because Pj — 4>{-^^) on by (|23l) . (Note 
that (f) equals 1 near 0, and so has an even extension to a function.) 
Finally, part (b) follows from Corollary [3] (b), once we observe that 

\i>jk,sM{^)\ = ^k\AM{x,^,k.B-^)l A,fc ~ B-^=. (60) 

This completes the proof. 



For notational simplicity, we take B = 2 for the rest of this article; the results would easily generalize to 
general B. 

Using Lemma [71 we obtain the following estimates on the norms of the ipj). . 
Lemma 11 For 1 < p < oo, we have 

||V',MAi||,-2-'-(i-'^- (61) 

Proof Let us call the estimate \\ipjk sM \\p — C2-'^^^p^ the majorization for this value of p, and call the 
reverse inequality the minorization. 

First we do the cases p = 1,2, oo. For p = oo, the majorization (by C2^ ) follows at once from (|60l) and 
Lemmalll For p = 1, the majorization (by C2~^) is Corollary [10] (b) . For p — 2, the majorization (by a 
constant) follows at once from the tight frame property. 

For the remaining estimates, we adapt arguments from [1 and 3 . Fix c>0, 0<i^<l such that 
> c on the interval [v, 1]. Using the orthonormality of the spin spherical harmonics, one obtains the 
minorization for p = 2 from 

\\i'3k.sM\\l=T.^^>'\Y^-(^3k)\'b'i^)>c2-'^ E ^ ^'>0- 

Im !'^<eis/43<l} 

The minorizations for p — l,oo now follow at once from the simple general inequality 

Il/ll2< ll/llill.flloo (62) 

and the majorizations for p ~ 1, oo. 

Thus we may assume 1 < p < oo, p ^ 2. The majorization follows from the general inequality 

11/11^ <ll/lli||/IIL"' (63) 
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and the majorizations for 1 and oo. 

For the minorization, we note that if g < 2 < r, and if < 6' < 1 is the number with 1/2 = 9/q+{l~6)/r, 
then one has the general inequahty 

11/112 < wfKwm-'. (64) 

If p > 2, the minorization fohows, after a brief computation, from (|64|) in the case q = 1, r = p, and the 
minorizations for 2 and 1. If p < 2, the minorization fohows, after a briefer computation, from ()64|) in 
the case q — p, r = 00, and the minorizations for 2 and 00. This completes the proof. 

5 Spin Besov Spaces and their Characterization 

The purpose of this section is the characterization of functional spaces by mixed needlets. This issue was 
already addressed by ( Ij, where the characterizations of Besov spaces by the asymptotic behaviour of 
spin needlet coefficients is addressed; here we aim at an analogous goal by focussing on mixed needlct 
coefficients. Most of our notations and of the arguments to follow are classical and close to those provided 
by ([!])■ We start by recalling that (compare ^T5\\ ) if 2^ > \s\, 

2^ 

^23, s = 'His , 

l = \s\ 

the space of spin functions spanned by spin spherical harmonics of degree up to 2^ . We let 

a,(i^,;p) inf ||F,-G.||^. , 

the error from the best approximations in that same space. The definition of Besov spaces is then natural. 

Definition 12 (^ll) (Spin Besov space) We say that the spin function G belongs to the Besov 
space of order {p, q, r; s} (written Fg € -B^'s) */ '^"■'^ only if 

a,{Fs;p) ^ e,2-^\ 
where {Sj} £ and p > 1, q, r > 0, s G N. The associated norm is 

\\Fs\\bp^, II-^sIIlp + IkjIU ■ 

Remark Recall that Lj is the largest integer with e^.s < 2^(-'+^\ and that by (1501) . Lj ~ 2^ as j — ?• 00. 
In particular, there is a c > such that 2^"'= < Lj < 2-'+'= for all j. Thus, if 2^^"" > \s\, 

Thus if we set 

a,(F,;p):= inf \\Fs - G^h. , 

it is evident that we could use the aj in place of the Cj in Definition [T^] to define the same spaces with 
an equivalent norm. 

The following characterization is provided by (^Ij) and extends to spin fiber bundles classical results on 
approximation spaces. 

Theorem 13 (JTl) If Fg G Lf, the following conditions are equivalent to Fg G BP^^: 
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1. 



2. 



3. 



jk,s\ \\^jk,s\\i^P 



IP 



£3j 



where {eij} , {e2j} , {esj} G i'', and Cp, Cp > 0. 

The fact that (1) implies (2) is easy, while the converse is standard from Hardy's inequahty. The fact 
that Fs € imphes (1) follows at once from using the dj in Definition [121 while the converse follows 

from Corollary [TU] (a) for the Pj, once one notes that for any Gs G Vlj_is, PjGs — 4>{-^tt)Gs — Gs, so 

\\P,Fs-Fs\\lp < \\PMFs-Gs\\l- + \\G,-F,\\lp^. 
The equivalence of (2) with (3) is established by first showing that for any F ^ L^, 

Cp\\Q,F,\\^.^ < jEl/^A/lkj-MllLsj <Cp{||Q,_iF,||^. + ||Q,J^,||^p + ||Q,+i^^,||^.} (65) 

We will prove a "mixed" analogue of (|65l) in Theorem [15] below, by a proof which is very close to the 
proof in HQ . 

Remark 14 In the previous Theorem, the crucial result is of course provided by (3), which provides the 
characterizations of Besov classes by means of the decay of spin needlet coefficients. This feature could 
be provided by many alternative formulations; in particular, as in the mixed case, one has (Ifl) that 



2^(1-1). 



The previous result can hence be formulated as follows: The measurable spin function Fs belongs to the 
Besov space of order {p, q, r; s} if and only if 



§2 



\Fsix)\Pdx 



i/p 



q/p' 



1/9 



< oo 



Our main result in this section is to show that the mixed needlet coefficients can play exactly the same role 
as the spin coefficients in the characterization of functional spaces, despite their different mathematical 
features. More precisely, we have the following alternative characterization of Besov spaces: 



Theorem 15 The function Fs G if belongs to the spin Besov space B^?^ if and only if 

i/p 

> =£4,2-^'-. 



IP II / IIP 



-^jk,sM\ \\^jk,sM\\ijP j 

where {s^j} G and Cp,Cp > 0. Equivalently, Fs £ B^9s */ '^'^'^ only if 



(66) 



\Fsix)\''dx 



1 i/p 



+ 



q/p' 



jk,sM I 



1/9 



< OO 
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Proof. Given Theorem [T^ and the results estabhshed in the previous section, the proof is rather standard 
and very close, for instance, to the arguments in ([T). 
By Theorem [131 it suffices to establish that for all F € if , 



(67) 



In addition to Lemma [TTl we will need the inequality 



sM (^) I < C'm E ■ 



Cm2^ 



1 M 



V{l + 2^d(a;,e,,)} 
which follows from the properties of e-nets (see [3] or [T]). 

To establish the rightmost inequality of (|67|. we note first that, in view of the mixed analogue of (|45)) . 

fc fc 

(In fact, the sums are clearly termwise equal at least for Fg G C^; the equality for general Fg e 
follows by a density argument.) 

The result will then follow if we can prove that, for all Gg E 

In view of Lemma [TT] and ([BS]) , the result can be established along exactly the same lines as in Lemma 
14 of (P); more precisely, by Holder's inequality 



§2 



1/p 



1-1/p 



whence 



k ° V"'^' k J 



IIP-1 II , IIP 

jk,sM\\iji \\Tjk,sM\\ i^p 



as desired, by Lemma [TT] and ([68[) . The proof of the rightmost inequality is hence completed. 

As far as the leftmost inequality of (|67| is concerned, again our arguments are very close to [Ij, Lemmas 

15 and 16. Note first that, by (|4ip and Lemma [U 



jk,sM ' 



at least for Fg E L^. Now using Holder's inequality 



E (-^s, V'jfc.sAl) ^jk,sM 
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- / {y2\(^'''^Jk,sM)\\i'jk,sMi^) 
•^S' I k 

- J^,^^\{^s^^3k,sM)\''\i^jk,sMi^)\^^^\'^jk,sMix)\^ dx 

< CV(P-'^Y.\(P^''^^k.sM)^k,k,sMh^<CY.\{Fs,^^^^^ , 

k ° k 

again by ((68|) and Lemma ITT] This gives the leftmost inequahty of (j67|) for i^; G L^, and hence for 
Fs £ LP if p > 2. If instead p < 2, and Fg G if is general, we take a sequence F"^ € approaching i^g 
in Lf , consider that inequality with F^ in place of Fs, and take the limsup of both sides in m, to obtain 
the inequality for Fs, as desired. ■ 

Theorem [TSl could be formulated more directly as: The section Fg belongs to the spin Besov space B^^g 
if and only if there exist {sj } G i'^ such that 

Ei/3.MA.r=e,2-^+^(^-i)}. 

k 

Combining Theorems [15] and I13[ we have the bounds 
where the i'^ sequence {e2j} is such that 

More explicitly, the asymptotic behaviour of the norms of needlet coefficients is of the same order for the 
spin and mixed spin case, despite the fact that the coefRcients in the two cases have a rather different 
nature (the {/3jfc^s} spin- valued, while {f3jk,sM} complex valued scalars). An alternative way to 
formulate this conclusion is the following. Define B^'^ (C) as the Besov space of complex- valued functions 
on the sphere. Then: 

The spin function Fg belongs to the spin Besov space BP?^ if and only if the scalar complex-valued 

function / = {/e + ifn) belongs to i?f'?(C) , i.e. 

{Fg e B^l) ^ {{fE + ifM) e i?r(C)} . 

Note that the complex- valued function / does not correspond to the function g introduced in (I2ip , indeed 
for a given array of coefhcients {aj^.s = + ''''^i^n\im ^® '^^^ write 



Ini Im V 



whence 



9{x) = E 4, ^l\ "-lm.sYlmix) ^ f{x) = y^^aim.sYlmjx) 

Irn V Im 



6 Statistical applications 

6.1 Estimation of angular power spectra and cross-spectra 

A major asset explaining the success of needlets for the analysis of cosmological data refers to their 
uncorrelation properties. More precisely, it was shown in ( 2 ), that for isotropic random fields, needlet 
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coefficients are asymptotically uncorrelated at any fixed angular distance as the frequency j diverges. 
This result was extended to the Mexican needlet case by [2I],|3B], and motivated many applications 
to astrophysical data, for instance (cross-)angular power spectrum estimation ([43], [11]), detection of 
asymmetries ([H]), bispectrum estimation ([30], [46], [45]) and many others. An analogous property was 
established for spin needlets in [1^; statistical techniques were then developed in 13 , while applications 
to CMB polarization data were detailed in [T2| . In this Section, we shall show how mixed needlets allow 
for further applications which have great physical interest and are not feasible by the pure spin approach, 
such as, for instance, estimation of cross-spectra between scalar and spin fields. 

To this aim, we shall focus on zero-mean, isotropic spin Gaussian random fields. As discussed by [14], [13], 
[32] and [33], the latter can be characterized by assuming that {a^,ahn} are complex- valued Gaussian 
random sequences satisfying 



and 



_™ , — ^1 



For m = 0, {a^jafcf} are real-valued Gaussian with the same moments. In the cosmological literature, 
the sequences {C;^, C*^} are known as the angular power spectra of the E and M modes; clearly in the 
Gaussian case they encode the full information on the dependence structure of the random field. In these 
area of applications, data are collected also on a standard scalar field (the so-called temperature of CMB 
radiation), which is again assumed to be Gaussian and isotropic with angular power spectrum 

p„T -T ^ s:l' s:m' r<T J7„T -E _ a'' A™' r-^-E p T -A/ _r/'cm'^TAf 

The cross-spectra {Cf^, Cf^^} are themselves of great physical relevance. The former is used to constrain 
cosmological parameters, in particular the so-called reionization epoch, while a detection of a non-zero 
value for the latter would entail an (unexpected) violation of parity invariance at the cosmological scales. 
Note that, for jointly isotropic random fields T, E, we have 

Cf^ = iJaL^L = E { [Re(aL) + » M^D] [Re(af„) - t Im(af„)] } 

= E {Re(aL) Re(af„)} + E {im(aL) Ini(aL)} + {Re(aL) Ini(aL)} - ^E {im(aL) M^Fm)} 
= E {Rc{afj Re(aL)} + E {im(aL) Ini(aL)} = , 

because {Re(a^), Im(a^)} = {im(a^), Re(a^)} by isotropy, where = denotes equality in distribution; 
the latter property follows (as in [5] and [37j, compare Theorem 7.2 in [T3]) from 

where {-DLim2 (-^)}mi denotes as usual the family of irreducible unitary representations of 5*0(3) by 
means of Wigner's matrices ([SD], [IS])- It follows, in particular, that the cross-power spectrum is always 
real-valued. 

We shall now provide an uncorrelation result that generalizes [2] and [13] . 

Theorem 16 Let Eg and T be jointly isotropic spin and scalar random fields (respectively), with angular 
power spectra such that 

Of = 5b(0'""", ct' = gA/(Or"-, Of = 5t(0^""^ 



aE,aM,aT > 2 , 



< CiU \ Ci > , i = 0, 1, 2, ... 



Assume also that for sufficiently large I 



\gE{l)\AgM{l)\,\9T{l)\>c>0 
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Then for all t > there exist > such that 

\Corr{l3j^^^.j^,l3^^^^.j^)\ , \Corr{(3jf^^.j^,j, l3^^^^.j^j)\ , \Corr{l3 ^^^.j., < 



Cr 



{l + 2MU,e,fcJ}' 



Proof. For the coeiEcents of the scalar random field T, the proof was provided in [2] (see [M] for the 
extension to the spin case). In view of the expressions provided in Section 4 for the needlet coefficients 
jk;E 1 1^ jk;M) ^'^'^ ^hc discussion following equation ([FT]), the proof is identical to the argument for the 
scalar case, and it is hence omitted for brevity's sake. ■ 

As we mentioned, mixed needlets allow for statistical procedures which were unfeasible in the scalar and 
pure spin cases. Consider for instance the issue of testing for a non-zero value of the cross-spectrum Cf^ , 
which is one of the main objectives of the ESA satellite mission Planck. Let us introduce the estimators 



We have easily, for A = E, M 



EfJ^ 



I. hmi l2m2 



)b(- 



/ei2s 



)aiimi-Tai2m2;Ey^^yhmi{^-jk)Yl2m2{£,jk)^jk 



Likewise 



Yar 



{fr} = {i:'*(# 



/iiI, 2Z + l 
47r 



cjct 



B3 ' An 



The following result is straightforward: 
Lemma 17 As j oo, for A — E,M 



iV(0,l) 



Proof. It suffices to note that 



j 



Bi 



and the summands satisfy all the assumptions of the classical Lindeberg-Levy Central Limit Theorem, 
see [11] for an analogous argument. ■ 

Remark 18 As in 'A3j , it is indeed possible to prove stronger results than Lemma \r^ namely, it can be 
shown that the same limiting result holds, if the estimator is constructed by using only the coefficients 
belonging to a connected subset ofS^. This result is important for astrophysical applications, where obser- 
vations are available only on subsets of the sphere, due to various forms of astrophysical contamination. 

Remark 19 It would be straightforward to exploit the properties of mixed needlet coefficients for many 
other statistical applications. To provide an example, it is possible to advocate estimation of the joint 
bispectrum of scalar and spin random fields, along the lines of the procedures advocated for the scalar 
case by fSO^ (see also \34^ , JSSjj, JSTJ). While these extensions are straightforward from the mathematical 
point of view, and hence omitted here for brevity 's sake, they are certainly of great practical importance 
for applications to CMB datasets, as we discussed in the Introduction. 
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6.2 Spin Nonparametric Regression 

While the uncorrelation properties of needlet coefBcients have aheady been widely exploited in statistical 
inference, the characterization provided for spin Besov spaces entails even richer statistical opportunities 
which are still almost completely open for research (see also [23], [25]) for classical results on adaptive non- 
parametric regression, [29l [27| for optimal spherical deconvolution methods, [4] and [26] for some results 
on needlet-based shrinkage estimation for densities on the sphere, and |28| for adaptive nonparametric 
regression on vector bundles). 

We envisage, in particular, applications to spin nonparametric regression by means of mixed-needlets 
shrinkage, in the following sense. Consider the regression model 

= F,(X,) + £. , i = l,2,...,n , 

where {Xi}^^^ ^ are (deterministic or stochastic) locations on the sphere Xi G S^, Fg G 3^?^ is a 
deterministic section of the spin fiber bundle and ^ sequence of random observational errors, 

themselves spin s variables. From the point of view of applications, we have in mind measurements of so- 
called weak gravitational lensing effects (see for instance |6]); here, Fg is the shear induced by gravitational 
effects on the image of distant galaxies, and {sk} are observational errors, due for instance to the intrinsic 
variability in the shape of the galaxies. The aim is to reconstruct Fg upon observations on {Ys{Xi)}; this 
is the object of a number of ongoing challenges, detailed for instance in [S]. Assume {Xi}^^^ ^ make 
up an (approximate) sequence of cubature points; we suggest to estimate F^ by means of the shrinkage 
procedure 

Fs{x) := ^'^Pjk g^Tpjk.sMi^) ' 

jk 



where 
and 



Pjk,sM = l^jk,sMK Pjk,sM > ctn) , i„ ^- as n ^ oo and c> 
'Pjk.sM — — ^Ys{Xi)xl: jj^^^^[Xi) ~ ( Fs{x)ipjk,sMi^)dx ■ 



The Besov space characterization opens the possibility of investigating optimality properties (in the 
minimax sense) of the shrinkage estimator Fg over Besov balls BP?g{Q), where Q < oo, defined as 

Fs '■ \\Fs\\ BfJ^ < Q ■ 

A full investigation of these issues will be reported elsewhere. 
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